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A theory of gravity with a generic action functional and minimally coupled to N matter fields has 
a nontrivial fixed point in the leading large N approximation. At this fixed point, the cosmological 
constant and Newton's constant are nonzero and UV relevant; the curvature squared terms are 
asymptotically free with marginal behaviour; all higher order terms are irrelevant and can be set to 
zero by a suitable choice of cutoff function. 
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There is growing evidence that a generalized version of 
General Relativity, possibly including terms of higher or- 
der in curvature, may be asymptotically safe in the sense 
of [1]. This property hinges on the existence of a non- 
trivial Fixed Point (FP) of the Renormalization Group 
(RG) flow of the gravitational couplings. Using Exact 
Renormalization Group Equation (ERGE) methods, it 
has been shown that a FP exists in 4 dimensions within 
the Einstcin-Hilbcrt truncation (where only the cosmo- 
logical constant and Newton's constant are retained) and 
there are various arguments suggesting that this FP is not 
a mere artifact of the truncation [2,3]. Independent evi- 
dence for the existence of a FP comes from recent results 
in causal dynamical triangulations [4]. Since in the real 
world gravity does not exist in isolation, but rather inter- 
acts with many matter fields, it is important to establish 
that matter interactions do not spoil the FP. (In fact, for 
applications to physics it is not logically required that 
pure gravity exists as a quantum theory: it is only neces- 
sary that a theory of gravity and realistic matter exists). 
In this direction, positive evidence has come from the 
analysis of arbitrary matter fields minimally coupled to 
gravity in the Einstcin-Hilbert truncation [5], and from 
the study of the coupled scalar-gravity theory with ar- 
bitrary couplings of the form (f> 2n and (f> 2n R with n > 
(thus including the original Einstein-Hubert terms) [6]. 

One of the central issues of this approach is the ex- 
tension of the results to include higher gravitational cou- 
plings. The beta functions of four-derivative gravity have 
been the subject of several papers in the past [7], also [8] 
and references therein; for the state of the art on this 
subject see [9]. Analyzing the structure of divergences in 
dimensional regularization, it has been shown that the 
couplings of the terms quadratic in curvature are asymp- 
totically free; however, due to the fact that dimensional 
regularization is not well-suited to the discussion of mass 
and threshold effects, the behaviour of the cosmological 
constant and Newton's constant has not been thoroughly 
understood. In [10] a term proportional to R 2 was in- 
cluded in the truncation of the ERGE and it was shown 



that its presence does not spoil the existence of the FP. 
In fact, the FP values of the cosmological constant and 
Newton's constant are only slightly changed relative to 
the Einstein-Hubert truncation, and the FP-value of the 
new term is quite small. 

In this note we shall consider a different approximation 
scheme, namely the large N expansion, where A" is the 
number of matter fields. This approximation has been 
first applied to gravity in [11] and was later used in [12] 
to prove the existence of a FP. I will use the ERGE in 
conjunction with the large A" approximation and the heat 
kernel expansion to rederive and extend this result. The 
truncation will consist of assuming that all matter self- 
interactions can be neglected. Nothing will be assumed 
a priori about the gravitational action. In the leading 
order in A" the analysis is so simplified, that we will be 
able to prove the existence of the FP for arbitrary grav- 
itational couplings, and furthermore we shall prove that 
there exists a scheme {i.e. a choice of cutoff function) 
such that all the coefficients of terms of cubic and higher 
order in the curvature vanish. 

The ERGE can be written, slightly symbolically, 
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where is the coarse-grained effective action, $ are 
all the fields in the theory, Tr denotes a trace over all 
degrees of freedom of the theory (sum over indices and 
integration over momenta), Rk is a cutoff function that 
suppresses the propagation of modes with momenta lower 
than k and t = log A:. The RG flow of the theory can be 
computed by making an Ansatz for the effective action 
Tfe, inserting it in |T]) and reading off the beta functions 
of the couplings. 

For simplicity we will begin by considering gravity cou- 
pled to N real scalar fields. In the spirit of an effective 
field theory, we will consider arbitrary gravitational ac- 
tions of the form 
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where O^ 1 ' = J d n x yfgM^ and M.\ n ' are polynomials 
in the curvature and its derivatives involving 2n deriva- 
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tives of the metric. The index i is used to label dif- 
ferent polynomials having the same n (and hence the 
same mass dimension 2n) . The first two polynomials are 
simply = 1 and J\4^ = R. The corresponding 

couplings are g^ = 2Z g k, where A is the cosmological 
constant and <jrW = —Z g = — 16 * G , where G is Newton's 
constant. For n = 2 I will follow [7-9] and choose as 
independent polynomials 
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where C is the Weyl tensor and G is 32ir 2 times the in- 
tegrand of the Gauss-Bonnet topological invariant. The 

(2) (2) 

coefficients g{ and <7g are sometimes denoted 1/A and 
— w/3A [8,9]. The second and fourth terms are total 
derivatives and do not contribute to local physical pro- 
cesses. We will assume that the scalar fields are all mass- 
less and minimally coupled to gravity, with inverse prop- 
agator z = —V 2 . For the running effective action Tk we 
assume that it is the sum of the scalar action and 

(n) 

where all couplings g\ are supposed to be fc-dependent. 
This defines our truncation. We will discuss later the sta- 
bility of this truncation against quantum fluctuations. 

In the 1/N expansion one assumes that the number 
of matter fields N — ► oo [11]. Here we shall restrict our- 
selves to the leading order approximation, where one sim- 
ply neglects the gravitational contribution with respect 
to the matter contribution. This eliminates most of the 
complications that arise in the gravitational ERGE cal- 
culations and provides the simplest way of arriving at a 
gravitational FP. Taking into account the cutoff function, 
the modified inverse propagator is Pk(z) — z + Rk(z); we 
can then write the ERGE (TTJ 
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The trace of a function / of the Laplacian can be writ- 
ten as 

/>oo 

Tr/ = V/(A l )=/ dtTrK(t)f(t) (5) 
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where Xi are the eigenvalues of the Laplacian, / is the 
Laplace transform of / and TvK(t) = J2 n e~ tx " is the 
trace of the heat kernel of —V 2 . (We assume that there 
are no negative and zero eigenvalues; if present, these will 
have to be dealt with separately.) 

One can then use the known properties of the heat 
kernel to calculate the trace in various circumstances. If 
we are interested in the local behaviour of the theory (i.e. 
the behaviour at scales k much smaller than the typical 
curvature) we can use the asymptotic expansion 

TrK(t) « Bat' 2 + B 2 r x + B 4 +B 6 t + ... (6) 
where B 2 k = J d A x y/g b 2 k, b 2 k being well-polynomials in 
the curvature tensor and its covariant derivatives. Then 
we get 
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We can rewrite these integrals in terms of the original 
function /: for n > 1 it is a Mcllin transform 
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where denotes the n-th derivative of /. Using {Zj, 



d t T k 



N 1 

y (4tt) 5 



d 4 Xy/g 



d t Pk 
Pk. 



\RQi 
6 



+ 



1 

180 V 2 



-C 2 



1 



-G 



d t Pk 
Pk. 



-R 2 + 6V 2 i? Q 



d t Pk 
Pk 



(11) 



+ 



Define dtg!- = pf 1 = fc 4 ~ 2 ™a,[™' ) . The coefficient of each where the numbers a\ can be read off (fTTj) . For all 

term in (jllj) is the beta function of the corresponding n^2 this leads to a fixed point 
coupling. The beta functions of the dimcnsionless cou- 
plings g^ 1 = gy^'k 271 ^ 4 are given by 

9tt ] = = (2- - 4)#> + aP (12) ^ = J^fP ^) 
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For n = 2 one gets instead the following solution: 

gf\k)=gf\k a ) + af ) Hk/h ) ) (14) 

As noted in [1], the condition of reaching a FP has to 
be imposed on those parameters for which a divergence 
would lead to divergent physical reaction rates. In the 
Wilsonian approach used here, it is enough to consider 
the tree level formulas deduced from the effective action. 
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In the case of the terms M\ and Ai\ z> , the interac- 
ts) (2) 

tions are not proportional a\ and , but rather to 
their inverses (in the notation of [8.9], they are given by 
A and — 3\/u>). The logarithmic running of the a^ 2 ^'s 
implies a FP. more specifically asymptotic freedom, for 
these coupling constants. The physical significance of the 
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parameters a 2 and a\' , which are the coefficients of to- 
tal derivative terms, is less clear in this respect, and we 
simply observe that the running derived above implies 
asymptotic freedom for their inverses. 

Choosing a cutoff function Rk determines the num- 

(n) 

bers a] and hence affects the position of the FP. This 
choice does not affect physically measurable quantities, 
however. One particularly convenient choice of cut- 
off function is the so-called optimized cutoff Rk{z) = 
[k 2 — z)0(k 2 — z), where 9 is the Heaviside step func- 
tion [13]. With this choice we get Q 2 (%^) = k 4 , 

Qi (%^) = 2k 2 , Q = 2 and Q n (%^) = 

for n < — 1. In terms of the cosmological constant and 
Newton's constant, the FP occurs for 
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The sign of Newton's constant is unphysical, but this is 
only a feature of the scalar fields (see eq.(18)). Perhaps 
the biggest surprise is the fact that in this scheme the 
FP-value of all the couplings with n > 3 is zero. This is 
due to the flatness of the function for small z. With 
the optimized cutoff, it would therefore be consistent to 
neglect all the terms with more than four derivatives of 
the metric. 

The dimension of the critical surface is determined 
by the linearized flow, which is defined by the matrix 



tJ — Since the in (| 1 2[) arc independent 

of the couplings, the eigenvalues of M are equal to the 
canonical dimensions. Therefore, the cosmological con- 
stant and Newton's constant are UV relevant (attrac- 

(2) 

tive), the g\ couplings are marginal and all the higher 
terms are irrelevant. 

With the same techniques, one can explore also other 
regimes. For large t, the heat kernel on a noncompact 
manifold has the asymptotic behaviour Kit) rts -nhvit~ 1 S 
where [15] 
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Since all the terms have the same dimension, they also 
have the same overall running. The beta functions of 
the operators appearing in S arc equal to (Aa Qit^), 
times the coefficient of the operator as it appears in 
S. Consequently, there is a FP where the FP-action 
is proportional to S, the exact coefficient being scheme- 
dependent. With the optimized cutoff, the FP-action is 
T* = 7nr-zk 2 S. We note that nonlocal actions of this 
form may have some relevance to cosmology [16]. 

Let us now consider the effect of other matter fields. 
Following [14], we can calculate the contribution to the 
ERGE of ns scalar fields, no Dirac fields, um Maxwell 
fields, all massless and minimally coupled. Unlike in [14], 
however, we shall not make any assumption about the 
background metric, so as to be able to discriminate the 
coefficients of various contractions of two curvatures. For 
each type of field we choose the cutoff function in such 
a way that the modified propagator has the form Pk(O), 
where C (s) = -V 2 on scalars, {D) = -V 2 + f on Dirac 
fields and {v) = -V 2 ^ + Ri* v on Maxwell fields in the 
gauge A = 1. We can then use the heat kernel coefficients 
given e.g. in [17]. In the case of the Maxwell fields, the 
calculation of the effective action with nonminimal oper- 
ators (corrisponding to other gauges) was given in [18], 
and shows that the results are gauge-independent. 

We assume that the numbers ns, nyy and tim are all 
of order N, or zero. The ERGE becomes 
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where the argument of each function appearing under the 
traces is the appropriate operator O and the last term is 
due to the ghosts. With the optimized cutoffs, the beta 
functions are characterized by the following coefficients: 
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We see that there is a FP for any value of ns, nw and 
tim- As expected, bosons and fermions contribute with 
opposite signs to the cosmological constant. Newton's 
constant is only positive if the ratio of scalar and Dirac 
to Maxwell fields is not too large. The couplings (<& ) -1 
are asymptotically free, tending to zero from above or 
below depending on the sign of a\ , and again = 
for n > 3. The sign of and a% guarantees that 
the Euclidean path integral is damped. Since this action 
represents the asymptotic UV behaviour of the theory, 
it is not required to reproduce the observed large-scale 
behaviour of gravity. 

Let us now discuss the consistency of the truncation 
that has been used in this paper. In general, one could 
expect processes with intermediate gravitons to gener- 
ate matter interactions even if they were initially set to 
zero. However, this will not happen in the UV limit if the 
interactions are asymptotically free. It has been shown 
in [6] that a large class of scalar interactions becomes 
asymptotically free in the presence of gravity. Thus, the 
truncation we used here for the scalar fields is probably 
consistent. Gauge fields are asymptotically free already 
in flat space and it is reasonable to expect that the situ- 
ation will not change in the presence of dynamical grav- 
ity. More troublesome are the Yukawa couplings, which 
in the standard model have positive beta functions. It is 
tempting to conjecture that the phenomenon that occurs 
for scalar fields is more general, and that there exists a 
"Gaussian-Matter FP" where all the matter interactions 
are asymptotically free and only the purely gravitational 
couplings have novanishing values at the FP. Earlier 
calculations reported in [8] suggest that this is the case. 
Then, the truncation adopted in this paper with regard 
to the matter fields would be a consistent truncation. 

The use of the large N limit has several advantages 
over previous calculations. First of all, it dramatically 
simplifies the derivation of the beta functions, making 
the result very trasparent. Furthermore, the calculations 
presented here are free of the gauge fixing ambiguities 
that occur when graviton loops are taken into account. 
This shows that the FP cannot be an artifact of the gauge 
choice. The most important result derived here is the 
existence of the FP for all the terms of higher order in 
curvature. What is perhaps even more striking, all the 
couplings with n > 3 can be made to vanish by a simple 
choice of cutoff scheme. It is known that the position 
of the FP is scheme dependent, but in the case of the 
cosmological constant and Newton's constant no sensible 
scheme is known were they vanish at the FP. 

One of the problems with applications of the ERGE is 
the absence of systematic ways of estimating the errors 
due to the truncation. If one can work with a sufficiently 
general consistent truncation, the large N limit intro- 
duces just such an expansion parameter. Calculating the 
next -to-leading order in l/N requires the evaluation of 
graviton contributions. This is obviously impossible with 



the full gravitational action ([2]), but the results presented 
here suggest that in certain schemes it may be consis- 
tent to set to zero all terms with n > 3. This would 
not be an a priori, uncontrollable truncation: due to the 
scheme-independence of the physical results, it would be 
equivalent to calculating with the full action. 
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